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The exactly solvable eigenproblems in Schrodinger quantum mechanics typ- 
ically involve the differential "shift operators" . In the standard supersym- 
metric (SUSY) case, the shift operator turns out to be of first order. In 
this work, I discuss a technique to generate exactly solvable eigenprob- 
lems by using second order shift operators. The links between this method 
and SUSY are analysed. As an example, we show the existence of a two- 
parametric family of exactly solvable Hamiltonians, which contains the 
Abraham-Moses potentials as a particular case. 



1. Introduction 

The number of exactly solvable eigenproblems in non-relativistic quantum mechanics is 
small, and most of them can be dealt with the factorization method. This technique, 
introduced long ago by Schrodinger,^ was analysed in depth by Infeld and Hull,^ who 
made an exhaustive classification of factorizable potentials. Later on, Witten noticed the 
possibility of arranging the Schrodinger's Hamiltonians into isospectral pairs (supersym- 
metric (SUSY) partners).^ The resulting supersymmetric quantum mechanics catalysed 
the study of hierarchies of 'exactly solvable Hamiltonians'. An additional step was Miel- 
nik's 'atypical' factorization^ through which the general SUSY partner for the oscillator 
was found; this technique was immediately applied to the hydrogen potential.^ Meanwhile, 
Nieto^ and Andrianov et. al7 put the method on its natural background discovering the 
links between SUSY, factorization and Darboux algorithm. These developments caused 
the renaissance of factorization and related algebraic methods, with particular attention 
focused on the first order differential shift operators. 

As can be noticed, however, the scheme is still narrow. An obvious generalization 
arises when higher order differential shift operators are used to connect the Hamiltonian 
pair. The idea of HSUSY (higher order SUSY), recently put forward by Andrianov et. 

17-18 ^ggg g^igQ ig-j^ incubated since 70-tieth.^''~^^ In this paper we will restrict ourselves 
to the case when the shift operator is of second order, and we name it SUSUSY. 
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2. Second order shift operator technique 

We postulate the existence of a second order differential operator interconnecting two 
different Hamiltonians if, H: 

HA"^ = A'^H, (1) 

Equality (1) imposes some restrictions to the functions {V{x),V{x), I3{x),^{x)y. 

V = V + 2(3', 2V + 6 = -2-f- (3', V" + (3V' = 2-f(3' - Y , (4) 

where 5 is an integration constant. We shall suppose that (3{x) is given and we shall 
express the other functions {V{x), V{x), 7(2;)} in terms of it. If we solve 7(0;) from 
second equation (4) and substitute the result in the third equation (4), we get: 

^ - 2/3'2 + (3' (3^ - (3(3" - 6(3' = (3V' + 2(3'V. (5) 
Multiplying by /3, it can be immediately integrated, yielding: 

''^'^ = 2^-[2p) -^+T + ^-2' 
where c is a new integration constant. The other two unknown functions become: 



Before going to the particular cases, let us notice a curious relation between the second 
order shift operator and Witten idea of the SUSY quantum mechanics. 

3. Second order SUSY (SUSUSY) 

According to Witten,'^ SUSY arises by defining a set of operators Qi that commute with 
the (supersymmetric) Hamiltonian Hg, 

[g„if,]=0, i^l---N, (8) 

and satisfy the algebra 

{Qi,Qj} = S,jHs, (9) 

where [-, ■] represents the commutator and {■, ■} the anticommutator. Now, with the aid of 
the operators A'^ , A of the previous section, one can construct a case of the supersymmetric 
algebra (8-9) with N = 2. 
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With this aim, define the supercharges:^'' 

where is given in (3) and A is its adjoint. Notice that Q'^ — = 0. Let us construct 
an operator, which we cannot abstain to call the SUSUSY 'Hamiltonian': 

Using the SUSY languaje, = A^A and H~ = AA^ should be called the SUSY partners. 
Notice that the SUSUSY 'Hamiltonian' Hgs commutes with the two supercharges Q and 
gt. The SUSY generators Qi = (Qt + Q)/^/2, Q2 = (Q^ - Q)/iV2, and H^s satisfy the 
supersymmetric algebra (8-9). 

Notice that the SUSY partners H~ are now the fourth order differential operators. 
It can be shown that commutes with H and commutes with H. Hence, can 
be a certain function of H and a function of H. Indeed: 

H-'^iH+^-J -c, H-^(h+^-J-c. (12) 
A physical Hamiltonian Hg can be defined (compare with the recent ideas of i'''-^^) , 

H.= ('l °), (13) 



H 

and the SUSUSY 'Hamiltonian' Hgs is related to Hs by means of: 

Hss= (^^ + (14) 

Thus, the SUSUSY 'Hamiltonian' Hg^ is a quadratic form of a physical Hamiltonian H^. 
The diagonal elements of Hg are the two Hamiltonians H, H of the previous section, 
which are related by the second order differential operators A, A'' (compare ^^~^^). 
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4. Example: the SUSUSY oscillator 

We shall now look for the SUSUSY analogue of the oscillator Hamiltonian 

We will try to show the existence of a 2-parametric family of potentials isospectral to 
V{x) = x^. This has to do with the general solution I3{x) of equation (6). This solution 
of course should include the ladder operator — {a^Y, where a) — —d/dx + x is the 
standard ladder operator of the harmonic oscillator. This means that for V{x) — x^ 
and I3p{x) = —2x. cqiiation (6) should become an identity, which fixes the constants to 
c — 1,S — 4. Substituting these results again in (6) and multiplying by 2/?^, we get: 

Pp" - ^ - 2/5^/3' + ^ - 4/3^ - 2x^(3^ + 2 = 0. (16) 

Let us notice the existence of an explicit solution more general than f3p = —2x. It 
arises after multiplying the standard raising operator by the operator of atypical 
factorizations,'' i.e., = b^a^ , leading to: 

The general solution of (16), which depends on two constants, should reduce itself to (17) 
as one of them takes a particular value (or one of them becomes a function of the other 
one). 

Here, I would like to present some partial numeric results obtained when (16) is inte- 
grated to provide /3(x) which arises as continuous deformations of the particular solutions 
(17). With this aim we choose the initial point (/3(0), /3'(0)) on the Poincarc plane close 
to (/3p(0), /9p(0)) and use then a standard numeric integration package^ to find i3{x) and 
to look for singularities in the corresponding SUSUSY potential V{x) = x^ + 2f3'{x). If 
there is no singularity, wc increase slightly /3'(0) maintaining /3(0) fixed, and repeat the 
integration until finding the upper threshold of /3'(0): above this threshold a singularity 
arises while below it disappears. A similar procedure is used to find the low threshold. 
After that, we make a small change of /3(0) along (/3j,(0), /^^(O)) and start again the whole 
process. In this way we can split the /9/9'-plane into the region where the SUSUSY poten- 
tial V{x) is free of singularities and the rest. Notice that the points (/3p(0), /3p(0)) provide 
a curve on (3 (3' -plane: 

/3;(0) = -2 + /3p2(0), 1/3,(0)1 <-^. (18) 



*We have employed the routine 'NDSolve' of 'Mathematica' 
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Departing from (18) we made the classification on figure 1. 
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Fig. 1. Classification of tlie /?/5'-plane into regions where tlie SUSUSY potential 
V{x) has no singularities (the shadowed regions) and those with singularity (the 
white regions). 

As we can see, there is a non-trivial region, the shadowed one, where the SUSUSY 
potential V{x) has no singularity; it comprises the curve in (18). For the purposes of this 
paper, to show that the general family of SUSUSY oscillator potentials is two parametric, 
our calculation brings already some information. 

The SUSUSY potentials V{x) for various values of the pair (/9(0), /9'(0)) lying at the 
shadowed region (no singularity, we have fixed /3(0) = —0.7) are shown in figure 2^. This 
family is richer than the Abraham-Moses (AM) SUSY potentials:"^'^^ 

This is so because V{x) is essentially two-parametric while (19) is just one-parametric. 
Indeed, the Vx{x) of (19) can be numerically reconstructed by solving (16) with the points 
of (18) as the initial conditions. The SUSUSY family obtained by this procedure coincides 
with a plot of the analytic results (19). 

'^Indeed, V{x) is displaced with respect to V{x) — a. quantity 5E = —4. This can be seen after 
reahzing that ^(a;) + 4 x^ when (/3(0), /3'(0)) (0, —2). Hence, we decided to represent on the vertical 
axis of figure 2 the potentials V{x) + 4 rather that V{x). 
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Fig. 2. The SUSUSY potentials +4 for some values of /3'(0) and /3(0) = -0.7. 
All pairs (/3(0), /3'(0)) lie in the region where there is no singularity for V{x). 

Interesting that the SUSUSY family V{x) embraces some cases of the widely discussed 
double well potentials (DWP). The dynamics of a system in these potentials is of some 
relevance, because it illustrates the differences between the classical and quantum regimes. 
In particular, it clearly shows one of the most intriguing quantum effects, the tunneling 
of the system from one well to the other as a result of the evolution. In most of the 
situations where a DWP is a SUSY pair of the oscillator potential, the DWP spectrum 
has one level more below the ground state energy of the oscillator. Moreover, it is usually 
symmetric with respect to a; = (see e.g. ^'^^). For our SUSUSY DWP, apparently, it is 
unneccessary to add any level below the ground state energy of the oscillator to generate 
the double well: the spectra of V{x) + 4 and V{x) = are equal. ^ The price to pay 
is that Vix) and Vix) = x^ are not precisely the SUSY partners, as shown in section 
3. Moreover, although V{x) is a double well, it turns out that it is not symmetric with 
respect to any point x = Xq. We hope that the SUSUSY treatment here presented can be 
implemented to other physically interesting potentials. 
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'^This is at the moment a hypothesis supported by our numerical plots of P{x) and V{x). (The continuity 
argument might be important.) 
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